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K*" Abstract 

\^ In interpreting the SNO experiments, accurate estimates of the ud reaction cross sections are 
O of great importance. We improve the previous estimates of our group by updating some of its 
^ ' inputs and by taking into account the resuhs of a recent effective-field-theoretical calculation. The 
CN ' new cross sections are slightly (~1 %) larger than the previously reported values. It is shown to 



be reasonable to assign 1% uncertainty to the pd cross sections reported here; this error estimate 
does not include radiative corrections, for which we refer to the literature. 
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1. Introduction 

The establishment of the Sudbury Neutrino Observatory (SNO) |T], ^ has motivated intensive 
theoretical effort to make reliable estimates of the neutrino-deuteron reaction cross sections [3- 
6]. One of the primary experiments at SNO is the measurement of the solar neutrino flux. By 
observing the charged- current (CC) reaction, Ved — > e~pp, one can determine the flux of the solar 
electron-neutrinos while, by monitoring the neutral-current (NC) reaction, v^d I'xPn {x = e, /i 
or r), one can determine the total flux of the solar neutrinos of any flavors. These features make 
SNO a unique facility for studying neutrino oscillations. SNO is also capable of monitoring the 
yield of the neutrino-electron elastic scattering (ES), u^e u^e, which also carries information 
on neutrino oscillations. The first report from SNO was concerned with the measurements 
of the CC and ES processes. By combining the SNO data on the CC reaction with the Super- 
Kamiokande data on ES strong evidence for z/g oscillations has been obtained 0. It is to 
be noted that the sharpness of this important conclusion depends on the precision of theoretical 
estimates for the i/rf-reaction cross sections. In the present communication we wish to describe 
our attempt to improve the existing estimates. 

^The SNO data on ES is consistent with the Super-K data but the latter has higher statistics. 
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We first give a brief survey of the theoretical estimates of the z/(i-reaction cross section that 
were used in the analysis in 0.0 A highly successful method for describing nuclear responses 
to electroweak probes is to consider one-body impulse approximation (lA) terms and two-body 
exchange- current (EXC) terms acting on non-relativistic nuclear wave functions, with the EXC 
contributions derived from one-boson exchange diagrams 0]. We refer to this method as the 
standard nuclear physics approach (SNPA)[^.0 The most elaborate calculation of the fd cross 
sections based on SNPA has been done by Nakamura et al. (NSGK) [^. Since the vd reactions in 
the solar neutrino energy {Ey < 20 MeV) is dominated by the contribution of the space component, 
A, of the axial current (A^), the theoretical precision of a^d is controlled essentially by the accuracy 
with which one can calculate the nuclear matrix element of A. Let us decompose A as A = 
AiA + Aexc, where Aja and Aexc are the lA and EXC contributions, respectively. Since Aja 
is well known, the theoretical uncertainty is confined to Aexc- Now, among the various terms 
contributing to Aexc; the A-excitation current (Aa) gives the most important contribution [ p!0[] , 
and Aa involves the coupling constants for the A^j_NA vertex, the vrA^A vertex and the pA^A 
vertex, and the corresponding form factors. Although the quark model is believed to provide 
reasonable estimates for these coupling constants, it is at present impossible to test their individual 
values empirically; only the overall strength of the A-excitation current can be monitored with 
electroweak processes in a few-nucleon system. NSGK therefore considered two methods for 
controlling the strength of the A-excitation current. In one method, by exploiting the fact that 
the A-excitation current features in the np^'yd amplitude as well, its strength is determined so 
as to reproduce the np 'yd cross section. The second method uses the tritium P decay rate, 
rf, and the strength of the A-excitation current is adjusted, as in Refs.[|10], [l^], to reproduce the 
well-known experimental value of rf . The first method was found to give a^d about 3% larger 
than the second method, and NSGK adopted this 3% difference as a measure of uncertainty in 
their calculation based on SNPA. 

Apart from SNPA, a new approach based on effective field theory (EFT) has been scoring 
great success in describing low-energy phenomena in few-nucleon systems [12-14]. Butler et al. 
(BCK) 10 applied EFT to the ud reactions, using the regularization scheme called the power 
divergence subtraction (PDS) [|15[. Their results agree with those of NSGK in the following sense. 
The EFT Lagrangian used by BCK involves one unknown low-energy constant (LEG), denoted 
by LiA, which represents the strength of A^-four-nucleon contact coupling. BGK adjusted Lia to 
optimize fit to the ai,d of NSGK and found that, after this optimization, the results of the EFT 
and SNPA calculations agree with each other within 1% over the entire solar-i/ energy region. 
Furthermore, the best-fit value of Lia turned out to be consistent with what one would expect 
from the ^^naturalness" argument [^]. The fact that the results of an ab initio calculation (modulo 
one free parameter) based on EFT are completely consistent with those of SNPA may be taken 
as evidence for the basic soundness of SNPA. 

Having given a brief survey of the existing theoretical estimates of ai^d, we now describe several 
points that need to be addressed for improving the estimates. We first note that, as pointed out 
by Beacom and Parke [|1^], the value of the axial coupling constant, qa, used in NSGK is not the 
most updated one. This obvious deficiency needs to be remedied. Secondly, in their treatment 
of Aexc; NSGK left out some sub-dominant diagrams, and therefore it is worthwhile to examine 



the consequences of using the full set of relevant Feynman diagrams |TT|. Furthermore, NSGK 



^ In what follows, a^^and cr^*^ stand for the total cross sections (in the laboratory frame) for the CC and 
NC reactions, respectively; in referring to a^^a,nd cr^'^ collectively, we use the generic symbol, a,yd- The incident 
neutrino energy in the lab-frame will be denoted by Ei,. 

^ This approach was called the phenomenological Lagrangian approach (PhLA) in H. 
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adopted as their standard run the case in which the strength of the A-excitation current was 
adjusted to reproduce the measured np jd rate. However, the np — >• d'y reaction governed 
by the vector current cannot be considered as a better constraint than rf for monitoring the 
effective strength of the A^NA vertex relevant to the axial-vector transition. In the present 
work, therefore, we adopt as our standard choice the case in which Aexc is controlled by rf . 
Thirdly, at the level of precision in question, radiative corrections become relevant [16-18]. In 
this communication, however, we do not address radiative corrections per se and simply refer 
to the literature on this issue |T^, |TB[. A related problem is what value should be used for the 
weak coupling constant. One possibility is to use the standard Fermi constant, Gp, which has 
been derived from yU-decay and hence does not contain any hadron-related radiative corrections. 
Another possibility is to employ an effective coupling constant (denoted by Gp) that includes the 
so-called inner radiative corrections for nuclear /3-decay. NSGK adopted the first choice. However, 
since the inner corrections are established reasonably well, it seems more natural to use Gp instead 
of Gp. We therefore adopt here Gp as the weak coupling constant (see below for more detail). 
An additional point that warrants a further study is the stability of the calculated value of 
against different choices of the NN interactions. NSGK investigated this aspect for a rather wide 
variety of modern high-quality NN interactions [O, EDI and found the stability of ai,^ at the 0.5% 



level. The interactions considered in NSGK, however, are all local potentials and have similar 
values of the deuteron D-state probability, Pd- Since the CD-Bonn potential [^] has a somewhat 
smaller value of Pd than the other modern high-quality A^A^ potentials, we study here whether 
the stability persists with the use of the CD-Bonn potential. 

Besides these improvements within the framework of SNPA, we present here a new comparison 
between SNPA and EFT. Park et al. [22-24] have developed an EFT approach wherein the 



electroweak transition operators are derived with a cut-off scheme EFT (a la Weinberg [12]) and 
the initial and final wave functions are obtained with the use of the high-quality phenomenological 
nuclear interactions. For convenience, we refer to this approach as EFT*. EFT* applied to the 
Gamow- Teller transitions contains one unknown LEG denoted by d^, which plays a role similar 
to LiA in BGK. In EFT*, however, one can determine d^i directly from Pf [^. This allows a 
parameter-free calculation of a,^di and very recently Ando et al. have carried out this type of 
calculation [ES. We present a comparison between our new results based on SNPA and those 



based on EFT*, and we argue that good agreement between them renders further support for the 
robustness of a^d obtained in SNPA. It will be seen that the new values of a^d are close to those 
given in NSGK, but that a significant improvement in error estimates has been achieved. 

2. Formalism 

We study the total and differential cross sections for the GG and NG reactions of neutrinos 
and antineutrinos with the deuteron: 

Ve + d ^e-+p + p (GG) (1) 

u^ + d ^v^ + n + p (NG) (2) 

Ve + d ->e+ + ?2 + n (z>-GG) (3) 

v^ + d ^v^ + n + p (z/-NG) , (4) 

where x = e, fi or r. We briefly describe our calculational framework and explain in what specific 
aspects we improve upon NSGK[pI. 
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The four-momenta of the participating particles are labeled as 



v/v{k) + d{P) ^ l{k') + N,{p',) + N^{p'^), (5) 

where ^ corresponds to for the CC reactions [Eqs.([l|),(§)], and to z/ or z/ for the NC reactions 
[Eqs.(0),(^)]. The energy-momentum conservation reads: k + P = k' + P' with P' = p[ + A 
momentum transferred from the lepton to the two-nucleon system is denoted by = k^ — k'^ = 
p'M _ ptJ^ the laboratory system, which we use throughout this work, we write 

F = {E,, k), k'^ = {E'„ k'), P^ = (M,, 0), P'^ = (P'°, P'), = {iv, q). (6) 

The interaction Hamiltonian for semileptonic weak processes is given by the product of the 
hadron current (Ja) and the lepton current (L^) as 



H'w' = '^Jdx[J^^ix)L^^'\x)+h.c.] (7) 

for the CC process and 



H^^ = ^Jdx[J^^{x)L^^\x) + \..c.] (8) 

for the NC process. Here is the weak coupling constant, and Vud is the K-M matrix element. 
For the weak coupling constant, instead of Gp = 1.16637 x 10^^ GeV~^ employed in NSGK, 
we adopt here Gp = 1.1803 x 10~^ GeV"^ obtained from 0+ 0+ nuclear /3-decays p6[.Q Gp 
subsumes the bulk of the inner radiative corrections.^ The K-M matrix element is taken to be Vud 
= 0.97400 instead of = 0.9749 used in NSGK. 

The leptonic currents, L^^^^ and L^*-^'^, are well known. The hadronic charged current is 
written as 

Jf^(=r) = Vt{x) + At{x\ (9) 

where V\ and Ax denote the vector and axial- vector currents, respectively. The superscript +(— ) 
denotes the isospin-raising (-lowering) operator for the z/(z/)-reaction. The hadronic neutral current 
is given by the standard model as 

J^^{x) = {l-2sm^9w)V^ + Al-2sm^9wV^, (10) 

where 6w is the Weinberg angle. is the isoscalar part of the vector current, and the superscript 
'3' denotes the third component of the isovector current. The hadron current consists of one- 
nucleon impulse approximation (lA) terms and two-body exchange current (EXC) terms. 

The lA currents are given in terms of the single-nucleon matrix elements of Jx- The standard 
parameterization for them is 

<Nip')\V^iO)\N{p)> = u{p')[fvlx + i^axpq']T^uip), (11) 
<N{p') I A±(0) I N{p)> = u{p')[fAlxl' + fpl'qx]T^u{p) , (12) 



The relation between Gp and the quantities used in is: Gp^ — (Gy/Kid)^(l + ^^), where is the 
nucleus-independent radiative correction. 

To be precise, the inner corrections for the CC and NC reactions may differ but the difference reported in the 
literature uM is comparable to the estimated uncertainty of our present calculation (see below). 
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where m is the average of the proton and neutron masses. For the third component of the isovector 
current, we simply replace with The isoscalar current is given as 

<Nip')\V^{0)\N{p)> = u{p')[fv7x + t^cTx,q']lu{p). (13) 

As for the dependence of the form factors we use the results of the latest analyses in 

Mql) = Gn{ql){l + ^^pV){l + Vr\ (14) 

fniql) = G'^(gJ)(/ip-/i„-l-/i„r7)(l + r7)"\ (15) 

fAiql) = -9AGA{ql), (16) 

27T7 

fliiQl) = GD{ql){^lp + ^ln-l + ^^n^l){'^ + ^l)~\ (18) 



with 



GD{qi) = 1- ..;; , (19) 




GAiql) = (1 1 , (20) 

^^''^ V 1.04GeVV 

where Hp = 2.793, /i„ = —1.913, rj = — and is the pion mass. For gA, we adopt the current 



standard value, (?a=1-267[29 , instead of (7a=1-254 used in NSGK. In addition, as the axial- vector 



mass in Eg. (pil|) , we use the value which was obtained in the latest analysis ||28|| of (anti)neutrino 
scattering and charged-pion electroproduction. The change in G^(g^) is in fact not consequential 
for aiyd in the solar-i/ energy region. Regarding fp, we assume PCAC and pion-pole dominance. 
A contribution from this term is known to be proportional to the lepton mass, which leads to very 
small contribution from the induced pseudoscalar term in our case. Although deviations from 



the naive pion-pole dominance of fp have been carefully studied|32], we can safely neglect those 
deviations here. For the lA current given above, we carry out the non-relativistic reduction in the 
same manner as in NSGK. 

We now consider the exchange currents (EXC). The axial- vector EXC, ^exc' consists of a 
pion-pole term and a non-pole term, ^exc- Using the PCAC hypothesis, however, we can express 
Aj^xc terms of the non-pole contribution alone: 

q^ - - 

^EXC = ^EXC 2 5"(9 ■ ^EXC - l^^EXCo)- (21) 

We therefore need only specify a model for the non-pole terms; the total contribution of ^exc 
be obtained with the use of Eg. (|21|) . Regarding the space component of the axial- vector current, 

we employ Aexc adjusted in such a manner that the experimental value of Ff be reproduced (see 
the discussion in Introduction). Following Schiavilla et al. 0], we consider the vr-pair current (de- 
noted by ttS*), p-pair current (pS*), 7r-exchange A-excitation current (Ayr), p-exchange A-excitation 
current (Ap) and vrp-exchange current (vrp). The explicit expressions of these two-body currents 
(acting on the i-th and j-th nucleons) are as follows. 

Aj(g;7rS) = %L {(^^ x r,)±^, x fc, - r± [q + z^, x (p, + p^)]} + (z ^ j), (22) 
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Ay {q-pS) = fA ml + k^- ^^'^''^ )xkj-i [a, x {a, x fc^- )] x (p, + ) } 

+ (tj X Xj)^ {qcr,-{crj x fc^-) + i(crj- x kj) x (p^+p-) - [cTj x (cTj x fc^)] x kj}) + {i ^ j), (23) 
_ , 1 « f2 rr ■ k 

AS(.;A.) = ^ (it - ^) ..^^ + fe| ' ^ ^ + - 

-± 4 gl{l + Kp)'^ Pn{kj) r + + 

A,j(q;Ap) = -^/a ^, . 2^1,2 {4^, (<T,xfc,)xfc, -(r, XT,)±cr,x[(^,xfe,)xfc,]} 

+ (25) 

Aj(q;^p) = 2/^^ ^_^Z_^2_^^/^(fe,)/4fc,)(r.xx,)± [(l + ^p)<T. x fe, - z(p,;+pO] + (* ^ j)-(26) 

Here mp and are the masses of the p-meson, and A-particle, respectively; /a is the axial form 
factor given in Eq.(|T6|). The total three- momentum transfer is g = fej + fcj, with fej^) being the 
momentum transferred to the i-th (j-th) nucleon; and p[ are the initial and final momenta of 
the 2-th nucleon. The form factors, f-jrik) and fp{k), for the pion-nucleon and p-nucleon vertices 
are parametrized as 

with = 4.8 fm~^ and Ap = 6.8 fm~^. The quark model has been used to relate the coupling 
constants of the vriVA, pNA and A^A^A vertices to the ttNN, pNN, and A^j^NN vertices, respec- 
tively. Schiavilla et al. have pointed out that the experimental value of rf can be reproduced 
if the strengths of A(A7r) in Eq. (p^) and A(Ap) in Eq.(p5D are reduced by a common factor of 
0.8. We employ here the same adjustment of A(A7r) and A(Ap). For the third component of the 
isovector current, we simply replace and (tj x Tj)^ with rf/2 and (tj x r^fjl, respectively. 
(The same prescription will be applied to the other exchange currents as well. ) For the time 
component we use the one-pion exchange current (the so-called KDR current ) , which gives 
the dominant EXC to ^o^jj- '^^^ explicit form of the KDR current, with a vertex form factor 
supplemented,^ reads 



Regarding the vector exchange currents, we first note that the time component should be 
negligibly small since its contribution vanishes in the static limit. For the space component, V, 
we take account of the pair, pionic, and isobar currents. As in NSGK, we adopt the one-pion 
exchange model for the pair and pionic currents and the one-pion and one-p-meson exchange 
model for the isobar current. Their explicit expressions are: 



y±(q;pazr) = -2^!y ( 1-Y f^ik. f''^^^' ^^\ t, x t,)^ + (^ ^ j), (29) 



Vt^iq■,p^on^c) = 2^ f ^) %.(fc.)/.(fc.) ^^^'A; (r. x r,)^ (30) 



^ For Aq and the vector currents, we use the same form factors as in NSGK. They are parametrized as in 
Eq.(p^), but the numerical values of and Ap are: A^r — 6.0 fm~^, Ap = 7.3 fm~^ 
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ir±f A\ -A fv + fM 



— jq X {cpkj X (cTj X kj)Tj + dpCTi x (fcj x (a-j x kj)){Ti x t^) | 



+ (31) 
The numerical values of the various parameters are 

^=0.08, com^ = 0.188, dim^ = -0.044, Cpm^ = 36.2, dp = -ic„. (32) 
47r 

As discussed in NSGK, these values lead to np — > dj cross sections that agree with the experi- 
mental values. 

Apart from the modifications explicitly mentioned above, the theoretical framework of the 
present calculation is the same as in NSGK ||5| and, for further details of the formalism, we refer 
to Sections II and III of NSGK 

3. Numerical results 

In reporting the numerical results of our calculation, we shall be primarily concerned with the 
standard case" ,f\ which is characterized by the following features. The calculational framework of 
the standard case is the same as in NSGK except for the specific points of improvements explained 
above. Also, the numerical values of the input parameters in the standard case are identical with 
those used in the standard run in NSGK, apart from the changes explicitly mentioned in the 
preceding section. As for the A^A^ interaction needed to generate the initial and final nuclear 
wave functions, the standard case employs the AV18 potential [|19|.0 In what follows, we largely 
concentrate on the standard case and discuss other cases (to be specified as needed) only in the 
context of assessing the model dependence. 

For the standard case, we have calculated the total cross sections and differential cross sec- 
tions for the four reactions Eq.(|l])-(^, up to Ei^ =170 MeV. In this communication, however, we 
concentrate on the quantities directly relevant to the SNO solar neutrino experiments and limit 
ourselves to the neutrino reactions (both CC and NC) for E^, < 20 MeV.^ The a^d corresponding 
to the standard case is shown in Table |I| as a function of The results given in Table |l| should 
supersede the corresponding results in NSGK. In the following we discuss comparison between the 
new and old estimates of a^a as well as error estimates for the new calculation. 

For clarity, when necessary, the total cross sections corresponding to the standard case of 
the present work are denoted by o"jyd{Netai), a^^(Neta.i) and cr^'"(Netai); those corresponding to the 
standard run in NSGK are denoted by a^dif^SGK), cr^'" (nsgk) and a^'"{NSGK). The ratio of (T^j"{Netai) 
to a^j" (NSGK) is given for several representative values of E^ in the first column of Table ^. Similar 
information for o"^*" is given in the second column. As the table indicates, (T^j^(Netai) is slightly 
larger than (T^j^(NSGK); the difference is ~ 1.3% for E^ 5 MeV, ~ 0.8% for E^ ~ 10 MeV, and 
~ 0.4% for El, ~ 20 MeV. A similar tendency is seen for as well. The origins of the difference 
between a^di^etai) and a^di^SGK) will be analyzed below. 



^ The standard case here should be distinguished from the standard run in NSGK. 

® The use of the AV18 potential here is consistent with the fact that, in |Q, the strength of d^ was determined 
using the 3-body wave functions obtained with the AV18 potential (and additional three-body forces). We also 
mention that the standard run in NSGK also uses the AV18 potential. 

^ A more extensive account of our calculation will be published elsewhere. The full presentation of the numerical 



results of the present work can be found at the web site: <http://nuc003.psc.sc.edu/~kubodera/NU-D-NSGK>. 

The limited precision of our computer code causes 0.1% uncertainty in a^d for < 20 MeV, apart from 
uncertainties due to the model dependence to be discussed later in the text. 
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Changing the weak couphng constant from Gp to Gp scales cjud by an overall factor of 
{Gp/GfY ~ 1.02. The effect of changing the value of qa can also be well simulated by an 
overall factor, since the vd reaction at low energies is dominated by the Gamow- Teller transition 
and hence a^d is essentially proportional to g\. Thus the change of qa from qa = 1.254 to qa = 
1.267 enhances a^d in the low-energy region by another factor of (1.267/1.254)^ ~ 1.02. 

In discussing the consequences of the change in Aexc^ it is convenient to introduce the terms, 
Models I and II. As described earlier, the standard case in the present calculation uses Aexc 
given in |TT| and recapitulated in the preceding section. We refer to this choice of Aexc as Model 



I. Meanwhile, Aexc used in the standard run in NSGK consists of A(A7r) and A(Ap) alone, 
and its strength is adjusted to reproduce the np 'yd rate. We refer to this choice of Aexc 
as Model II. For each of Models I and II, Table |^ gives the contributions from the individual 
terms in as well as that from V^, the vector current. The table indicates that in either case 
the corrections to the lA values are dominated by the contributions from A(A7r) and A(Ap). 
To facilitate further comparison between Models I and II, we consider the ratio, ^, defined by 
^ = [cr^rf(IA+AExc) ~ c";/d(IA)]/cT,^rf(IA). Here, a^dij-^) is the result obtained with the lA current 
alone, while (Ti,rf(IA+ Aexc) represents the result obtained with the lA current plus Aexc- Fig. |I] 
gives ^ for the CC reaction as a function of E^. The solid line shows ^ for Model I, while the dashed 
line gives ^ for Model II. It is seen that the contribution of Aexc in Model I is smaller than that 
in Model II by 2 ~ 4%. This difference is mainly due to the reduced strength of the A-excitation 
currents in Model I. For further discussion, we normalize ^(Modei ii) for the CC reaction by an 
overall multiplicative factor chosen in such a manner that the normalized ^ (Model ii) reproduces 
^(Modei I) at the reaction threshold. This normalized result is given by the dash-dotted line in 
Fig.|l[ We observe that the dash-dotted line exhibits a slight deviation from the solid line (Model 
I). This deviation reflects a slight difference in the i?,y-dependences of the Aexc contribution for 
Models I and II. The main cause of this difference can be traced as follows. From Table |^ we can 
deduce that the dominant contributions coming from the A-excitation currents have almost the 
same ii^j^- dependence for Models I and II, although their absolute values differ for the two models. 
On the other hand, the contributions of A(7r5'), A{pS) and A{7ip), which are included only in 
Model I, are virtually i?i/-independent, and their magnitudes are small. These features lead to the 
slighly weaker _E,^-dependence in AExc(Model I) than in AExc(Model II). The behavior of C, for 
the NC reaction (not shown) is similar to ^ for the CC reaction. 

The error estimate adopted in NSGK essentially consists in taking the difference between 
Models I and II as a typical measure of model dependence. As mentioned in the introduction, 
however. Model II, which fails to explain rf , should not be given the same status as Model I. 
To attain a more reasonable estimate of the theoretical uncertainty, we propose the following 
interpretation of the feature seen in Fig. |[ The fact that Model I has been adjusted to reproduce 
rf means that it can yield model- independent results at a specific kinematics but that, without 
additional experimental information, the ii^i,- dependence of a^d cannot be fully controlled. This 
uncertainty may be assessed from the difference between the solid and dash-dotted lines in Fig.^ 
From this argument we assign 0.2% uncertainty to the contribution of Aexc to a^d in the solar 
neutrino energy range, < 20 MeV. 

We recapitulate our discussion regarding the change from ^^.^(nsgk) to cryrf(Netai): a ~4% en- 
hancement of a^d due to the changes in the Fermi constant and qa and a ~3% reduction due to 
the use of the new Aexc (Model I) that reproduces rf . These two changes partially cancelling 
each other, the net result is the enhancement of (j,yd by ~1%, and this is what is seen in Table ^ 

As mentioned earlier, an additional important measure of reliability of our SNPA calculation is 
obtained by comparing it with the results of an EFT* calculation by Ando et al. p5l . By using the 
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value of the low-energy constant, d^, fixed to reproduce the experimental value of rf U^, Ando et 



al. 1 25 1 have carried out a parameter-free EFT- motivated calculation of a^d- Although the cut-off 



regularization method used in p5[ can introduce the cut-off dependence into the formalism, it has 
been checked I^Sf that this dependence is negligibly small for a physically reasonable range of the 
cut-off parameter; the relative variation in a^a only amounts to 0.02%, which is much smaller than 
the above-mentioned 0.2% uncertainty inherent in our SNPA calculation. In fact, the uncertainty 
in ai^d obtained by Ando et al. is dominated by the 0.5% error resulting from the uncertainty in 
the experimental value of rf . We now compare a^di^etai) with cr^di^FT*) obtained in the EFT* 
calculation of Ando et al. ||25[|. Since Ref.[^ only includes the s-wave of the final A^A^ state, we 
compare (Tyrf(EFT*) with a^d{s-waye) , which represents the s-wave contribution to a^d calculated 
for the standard case. The ratio, rj = ai,d{EFT*) / ai,d{s-wave) , is shown in Table 0, from which we 
can conclude that SNPA and EFT give identical results at the 1% level. 

We proceed to consider the A^A^ potential dependence. As mentioned, the CD-Bonn potential 
is somewhat distinct from the potentials considered in NSGK, in that it has a significantly smaller 
D-state probability; Pc){CD-Bonn) = 4.2% as compared with Pd(avi8) = 5.8%. We compare in Table 
|] the a^J^ obtained with the AV18 and CD-Bonn potentials. The difference between the two cases 
is found to be practically negligible. With the CD-Bonn potential, because of its larger S-state 
probability, the contribution from the lA current becomes larger, whereas the contribution from 
^Exc is smaller due to its reduced D-state probability. Our explicit calculation demonstrates that 
the cancellation between these two opposing tendencies is almost perfect, providing a yet another 
manifestation of the robustness of the calculated cXj^rf. A similar stabilizing mechanism was noticed 



by Schiavilla et al. in their study of the pp-fusion cross section [11 



Apart from the absolute values of a^^and (J^d i ^^e ratio, R = cr^d l^vd i is also an important 
quantity for SNO experiments. As emphasized by Bahcall and Lisi a measurement at SNO 
of the ratio of the number of the NC to CC events (the NC/CC ratio) would place stringent 
constraints on various neutrino oscillation scenarios. Since the precision of the predicted value 



of the NC/CC ratio is affected by the uncertainty in R (see Fig. 7 in |32|), we discuss the model 



dependence of R. Table ^ gives the values of R calculated for the various cases discussed above. 
The table indicates that the model dependence of R is smaller than that of a'^^ and them- 
selves. The simple lA calculation gives a value of R that agrees with Rstnd within 0.2%. The 
variance between the standard case and Model II is less than 0.5% - it is to be recalled that Model 
II is a rather extreme case. Furthermore, the difference between the standard case and the EFT* 
calculation does not exceed 0.4%. We therefore consider it reasonable to assign 0.5% accuracy to 
R. This is an improvement by a factor of ~2 over the precision reported in 0, 3. 

The stability of R can be understood as follows. We first note the following two features. (1) 
The contribution of the isoscalar current, which only participates in the NC reaction, is negligibly 
small in our case; (2) Although the iso-vector vector and axial-vector currents enter in different 
ways into the nuclear currents responsible for the CC or NC reactions, the contribution of the 
vector current is much smaller than that of the axial-vector current in the solar neutrino energy 
regime. As a consequence of these two facts, the transition operators for the NC and CC reactions 
in the present case are, to good accuracy, related by a rotation in isospin space. So, if there were 
no isospin breaking effects in the nuclear wave functions, the CC and NC transition amplitudes 
would be simply related by the Wigner-Eckart theorem in isospin space, leading to the complete 
model independence of R. In reality, there are isospin-breaking effects in the two-nucleon wave 
functions, but these '^external effects'^ are expected to be under good control so long as one uses 
high-quality A^A^ potentials that reproduce the A^A^ data accurately. 
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4. Discussion and Summary 



Although we do not directly address the issue of radiative corrections (RC) here, we make a few 



remarks on it. RC can affect CTj^rf at the level of a few percents. According to Kurylov et al. 
RC increases cr^j^ by 4% at low and by 3% at the higher end of the solar neutrino energy, while 
RC leads to an E'^-independent increase of by ~ 1.5%. The RC for a^^ consists of the "inner" 
and "outer" corrections. The former is sensitive to hadronic dynamics but energy- independent, 
while the latter is largely independent of hadronic dynamics but has energy-dependence. The 
use of experimental value oi Gp obtained from 0+ — *• O"*" nuclear /5-decays allows one to take 
account of the bulk of the "inner" corrections. To obtain reasonable up-to-date estimates of the 
remaining "outer" corrections, we may proceed as follows. For o"^'", we may adopt as the "outer" 
correction the difference between the result of Kurylov et al. (4% - 3%) and the estimated "inner" 
corrections (2.4%). For a^*", there is no "outer" corrections at the level of precision of this article. 
In adopting this prescription, we are leaving unaddressed a delicate issue of the possible difference 
between RC for the single nucleon and RC for multi-nucleon systems, but this seems to be the 
best we can do at present. 

To summarize, we have improved NSGK's calculation for the vd reactions by updating 
some of its inputs and with the use of the axial-vector exchange current the strength of which is 
controlled by Ff . We have also taken into account the results of a recent parameter-free EFT* 
calculation [^. The new value of aud, denoted by cri.rf{Netai), is slightly larger than (T,yrf(NSGK) 
reported in cr,^rf(Netai)/cri,d{NSGK) ~ 1.01. Based on the arguments presented above, we consider 
it reasonable to assign 1% uncertainty to crj^rf{Nctai) given in Table |l], and ~0.5% uncertainty to 
R. The results in Table |I|, however, do not include radiative corrections except for those already 
incorporated into the empirical value of Gp, which subsumes the bulk of the inner radiative 
corrections for nuclear /3-decay. With the inclusion of the remaining radiative corrections, a^j" is 
likely to become larger than (T^j"(Netai) by up to ~2 %, while a^^ is expected to lie within the 
quoted ~1 % error of (T^'"(Netai). 
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Figure 1: Contributions of Aexc to cr^^; C, defined in the text is plotted for Model 1 (solid line) 
and Model II (dashed line). The dash-dotted line represents the "normalized" version of Model II 
described in the text. 
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Table 1: Calculated values of a^^ and a^^ in units of cm^. The "-a;" in the parentheses means 
10-^; thus an entry like 4.579(-48) stands for 4.579 x 10"^^ cm^. 



ieY) 




e pp 




vpn 


E, 

(McV) 


Ved 


e pp 


vd 


vpn 


1.5 


4.680 


(-48) 


0.000 


(0) 


8.8 


1.911 


(-42) 


7.530 


(-43) 


1.6 


1.147 


(-46) 


0.000 


(0) 


9.0 


2.034 


(-42) 


8.070 


(-43) 


1.8 


1.147 


(-45) 


0.000 


(0) 


9.2 


2.160 


(-42) 


8.629 


(-43) 


2.0 


3.670 


(-45) 


0.000 


(0) 


9.4 


2.291 


(-42) 


9.209 


(-43) 


2.2 


7.973 


(-45) 


0.000 


(0) 


9.6 


2.425 


(-42) 


9.809 


(-43) 


2.4 


1.428 


(-44) 


4.346 


'-47) 


9.8 


2.565 


(-42) 


1.043 


(-42) 


2.6 


2.279 


(-44) 


4.322 


;-46) 


10.0 


2.708 


(-42) 


1.107 


(-42) 


2.8 


3.369 


(-44) 


1.478 


;-45) 


10.2 


2.856 


(-42) 


1.173 


(-42) 


3.0 


4.712 


(-44) 


3.402 


'-45) 


10.4 


3.007 


(-42) 


1.241 


(-42) 


3.2 


6.324 


(-44) 


6.372 


;-45) 


10.6 


3.164 


(-42) 


1.311 


(-42) 


3.4 


8.216 


(-44) 


1.052 


;-44) 


10.8 


3.324 


(-42) 


1.383 


(-42) 


3.6 


1.040 


(-43) 


1.594 


'-44) 


11.0 


3.489 


(-42) 


1.458 


(-42) 


3.8 


1.289 


(-43) 


2.274 


;-44) 


11.2 


3.658 


(-42) 


1.534 


(-42) 


4.0 


1.569 


(-43) 


3.098 


;-44) 


11.4 


3.832 


(-42) 


1.612 


(-42) 


4.2 


1.881 


(-43) 


4.072 


'-44) 


11.6 


4.010 


(-42) 


1.693 


(-42) 


4.4 


2.225 


(-43) 


5.202 


;-44) 


11.8 


4.192 


(-42) 


1.775 


(-42) 


4.6 


2.604 


(-43) 


6.492 


;-44) 


12.0 


4.379 


(-42) 


1.860 


(-42) 


4.8 


3.016 


(-43) 


7.947 


'-44) 


12.2 


4.570 


(-42) 


1.947 


(-42) 


5.0 


3.463 


(-43) 


9.570 


;-44) 


12.4 


4.766 


(-42) 


2.035 


(-42) 


5.2 


3.945 


(-43) 


1.136 


;-43) 


12.6 


4.966 


(-42) 


2.126 


(-42) 


5.4 


4.463 


(-43) 


1.333 


'-43) 


12.8 


5.171 


(-42) 


2.219 


(-42) 


5.6 


5.017 


(-43) 


1.548 


;-43) 


13.0 


5.380 


(-42) 


2.314 


(-42) 


5.8 


5.608 


(-43) 


1.780 


;-43) 


13.5 


5.923 


(-42) 


2.561 


(-42) 


6.0 


6.236 


(-43) 


2.031 


'-43) 


14.0 


6.495 


(-42) 


2.822 


(-42) 


6.2 


6.902 


(-43) 


2.300 


;-43) 


14.5 


7.095 


(-42) 


3.095 


(-42) 


6.4 


7.605 


(-43) 


2.587 


;-43) 


15.0 


7.724 


(-42) 


3.382 


(-42) 


6.6 


8.347 


(-43) 


2.894 


'-43) 


15.5 


8.383 


(-42) 


3.682 


(-42) 


6.8 


9.127 


(-43) 


3.219 


;-43) 


16.0 


9.071 


(-42) 


3.995 


(-42) 


7.0 


9.946 


(-43) 


3.562 


;-43) 


16.5 


9.789 


(-42) 


4.323 


(-42) 


7.2 


1.080 


(-42) 


3.925 


'-43) 


17.0 


1.054 


(-41) 


4.663 


(-42) 


7.4 


1.170 


(-42) 


4.308 


;-43) 


17.5 


1.131 


(-41) 


5.017 


(-42) 


7.6 


1.264 


(-42) 


4.709 


'-43) 


18.0 


1.212 


(-41) 


5.385 


(-42) 


7.8 


1.362 


(-42) 


5.130 


'-43) 


18.5 


1.296 


(-41) 


5.767 


(-42) 


8.0 


1.464 


(-42) 


5.571 


r43) 


19.0 


1.383 


(-41) 


6.162 


(-42) 


8.2 


1.569 


(-42) 


6.031 


;-43) 


19.5 


1.474 


(-41) 


6.571 


(-42) 


8.4 


1.679 


(-42) 


6.511 


'-43) 


20.0 


1.567 


(-41) 


6.994 


(-42) 


8.6 


1.793 


(-42) 


7.010 


;-43) 
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Table 2: Comparison of the present results with those of NSGK [^. The ratio, a^di^ctai) /a^di^SGK), 
is given for representative values of E^. 



El, (MeV) —f e pp vd —t vpn 

5 1.013 1.011 

10 1.008 1.006 

15 1.006 1.003 

20 1.004 1.001 



Table 3: For Models I and II are shown the cumulative contributions to cr^ from the various 
components in the current. The row labeled "lA" gives cr^j" obtained with the lA currents in 
and V^, and the next row labeled "+Vexc" gives cr^j" that includes the contributions of the lA 
currents and Vexc^ the exchange current in V . Similarly, an entry in the n-th row (counting from 
the row labeled "lA") includes the coherent contributions of all the currents listed in the first n 
rows. The numbers in the last row are obtained with the full currents. The parenthesized number 
in the n-th row gives the ratio, (T^*"(n-th row)/ a'^^{{n-l)-i}i row), which represents a factor by 
which cr^'" changes when the new term is added. 



^ud (xlQ-^^cni^) 



Model I 



E^ 5 MeV 10 MeV 15 MeV 20 MeV 

~TA 0.3397 ( - ) 2.646 ( - ) 7.526 ( - ) 15.23 ( - ) 

+Fexc 0.3401 (1.001) 2.654 (1.003) 7.560 (1.005) 15.33 (1.006) 

+7rA 0.3474 (1.022) 2.719 (1.025) 7.758 (1.026) 15.74 (1.027) 

+pA 0.3448 (0.992) 2.695 (0.991) 7.687 (0.991) 15.59 (0.991) 

+ttS 0.3456 (1.002) 2.702 (1.003) 7.707 (1.003) 15.63 (1.003) 

+pS 0.3447 (0.997) 2.694 (0.997) 7.682 (0.997) 15.58 (0.997) 

+TT-P 0.3463 (1.005) 2.708 (1.005) 7.724 (1.005) 15.67 (1.006) 

+Ag;xc 0-3463 (1.000) 2.708 (1.000) 7.724 (1.000) 15.67 (1.000) 



Model II 



E^ 5 MeV 10 MeV 15 MeV 20 MeV 

~TK 0.3397 ( - ) 2.646 ( - ) 7.526 ( - ) 15.23 ( - ) 
+Vexc 0.3401 (1.001) 2.654 (1.003) 7.560 (1.005) 15.33 (1.006) 
+7rA 0.3612 (1.062) 2.841 (1.071) 8.128 (1.075) 16.52 (1.078) 
+pA 0.3567 (0.988) 2.801 (0.986) 8.007 (0.985) 16.26 (0.985) 
+Ag,xc 0.3567 (1.000) 2.801 (1.000) 8.008 (1.000) 16.26 (1.000) 
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Table 4: Comparison of SNPA and EFT calculations. The ratio, 77 = (7,^^ (eft*) /(7i,(i(s— wave), is 
given for representative values of E^,. 



(MeV) 


Ved —^epp 


i^d — > upn 


5 


1.003 


1.004 


10 


1.001 


1.003 


15 


0.999 


1.002 


20 


0.998 


1.001 



Table 5: Dependence of cr^^ on NN potentials. 'Bonn' and 'AV18' represent the results obtained 
with the CD-Bonn and the AV18 potentials, respectively. 





5 MeV 


10 MeV 


15 MeV 


20 MeV 




Bonn AV18 


Bonn AVIS 


Bonn AV18 


Bonn AV18 


lA 


0.09459 0.09390 


1.091 1.083 


3.327 3.300 


6.871 6.814 


lA+EXC 


0.09557 0.09570 


1.104 1.107 


3.373 3.382 


6.973 6.994 



Table 6: The ratio, R = a{vd — > vnp) / a{v^d e pp), calculated for representative values of Ei,. 
The second column gives Rgtnd corresponding to the standard case. The third, fourth and fifth 
columns give Ria/ Rstnd, -R[Modei ii]/-Rsted and Rieft*]/ Rstnd, respectively. Here Ria corresponds to 
R obtained with the lA current alone, while i?[Modei n] and R[eft*] corresponds to Model II and 
EFT*, respectively. 



E„ (MeV) R.tnd lA Model II EFT* 

5 0.2764 1.000 1.004 1.001 

10 0.4087 1.001 1.004 1.002 

15 0.4378 1.002 1.004 1.004 

20 0.4464 1.002 1.005 1.004 
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